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Abstract 

The integrable open-boundary conditions for the Bariev model of three coupled 
one-dimensional XY spin chains are studied in the framework of the boundary 
quantum inverse scattering method. Three kinds of diagonal boundary K-matrices 
leading to nine classes of possible choices of boundary fields are found and the 
corresponding integrable boundary terms are presented explicitly. The boundary 
Hamiltonian is solved by using the coordinate Bethe ansatz technique and the Bethe 
ansatz equations are derived. 

Key words: Integrable spin chains; Algebraic Bethe ansatz; Yang-Baxter algebra; 
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1 Introduction 



Since tlie discovery of high temperature superconductivity in cuprates [1] , a 
tremendous effort has been made to uncover a theoretical framework capable of 
explaining this amazing phenomenon. It is a general belief that the properties 
of strongly correlated electron systems showing a non-Fermi liquid behaviour 
are closely related to those materials showing high Tc superconductivity. This 
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has caused an increasing interest in strongly correlated electron models [2- 
5]. Integrable models that have been widely studied include the ID Hubbard 
model solved by Lieb and Wu [6] and the supersymmetric t-J model [7,8]. 
Another model with special features relevant for high T — c superconductivity 
is the ID Bariev [9,10], as it exhibits the existence of hole pairs of Cooper 
type. 

On the other hand, important progress in the realm of completely integrable 
systems is the generalization of the usual quantum inverse scattering method 
(QISM) [11-13] to incorporate open boundary conditions that preserve in- 
tegrability [14-17]. The presence of the integrable boundaries leading to a 
pure back-scattering on each end of a quantum chain results in rich physi- 
cal phenomena [18-2 3] . It has been clarified that ID quantum systems with 
boundary fields are closely related to impurity problems. Also due to their 
connection to the Kondo problem and boundary conformal field theory in low- 
dimensional quantum many-body systems, the integrable boundaries could be 
a useful non-perturbative way to investigate the impurity effects in the con- 
densed matter physics. The prototypical ID Hubbard model and t-J model 
have been thoroughly investigated with boundary impurities and boundary 
fields [18-25]. Subsequently, other strongly correlated electron systems, such 
as the ID Bariev model of two coupled XY spin chains have been partially 
studied and the reflections equations (RE) solved[26]. In that case, a new 
class of boundary reflection K-matrices for the model lead to pure magnetic 
boundary flelds in the Hamiltonian which may have a feasible realization by 
applying boundary external fields in experiments on quantum wires. An in- 
teresting aspect here is that after Jordan- Wigner transformation, the coupled 
XY spin chains can be presented as correlated electron systems where the hop- 
ping terms depend on the occupation numbers at sublattices. This is thought 
to be useful in studying conductivity properties in such non-fermi liquids. 
Along this line, the integrability and boundary conditions for a Bariev model 
of three coupled XY chains has been studied in the frame work of the QISM 
by Zhou and coworkers [27,28] recently. However, we have good reasons to 
expect that the model would permit other kinds of integrable boundary terms 
associated with some new solutions to the RE. We revisit the model seeking 
a complete understanding of the open boundary conditions for the model. As 
noticed in [28] the quantum i?-matrix of the model we study does not possess 
the crossing-unitarity. This causes a violation of the isomorphism between K^- 
and /C_-matrix which satisfy two reflection equations (RE) separately. The in- 
tegrability at left and right boundaries demand rather complicated RE which 
in turn expose bulk symmetries on the left and right boundary terms. In order 
to maintain the integrability of models where the i?-matrix does not posses 
the crossing-unitarity property some new objects have to be introduced . It is 
found that the three coupled XY chains permits three kinds of K-matrices to 
each RE which lead to nine classes of integrable boundary terms containing 
the ones for the two coupled XY spin chain as special cases. These integrable 
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boundary terms containing different on-site Couiomb interactions and cliem- 
ical potentials reveal tlie symmetry of exclianging tlie sublattices. We derive 
tlie Hamiltonian of tlie model with nine classes of boundary terms from the 
expansion of the boundary transfer matrix around the zero spectral-parameter 
point up to the fourth, third and second orders respectively. Furthermore, we 
solve the boundary model by means of the coordinate Bethe ansatz method 
and derive the Bethe ansatz equations. 

The paper is organized as follows. In section 2 we construct the Bariev model of 
three coupled XY spin chains with nine classes of boundary fields by means of 
the QISM adapted to special boundary conditions. The basic quantities, e.g., 
the i?-matrix, the matrices K± defining the boundary terms, the monodromy 
matrices and the transfer matrices arc defined. In particular, we solve two RE 
separately and obtain three independent classes of solutions to each REs using 
the variable-separation method. The relation between the transfer matrices 
and the Bariev Hamiltonians with different boundary fields is established. 
Section 3 is devoted to the solution of the models through coordinate Bethe 
ansatz method. Section 4 presents our conclusions. 



2 Boundary K-matrices for the coupled spin chains 

We consider a spin chain model defined by the following Hamiltonian 

^=E^.,.+i + ^^^ + 4^ (2.1) 



where -Hjj+i denotes the bulk Hamiltonian density of three XY spin chains 
coupled to each other [27] 



^3,3+1 = E (^/(a)^j+l(a) + ^ 3{oc)^t+l{a)) (2-2) 
o 

exp[?7 ^ ^^+e{a'-a){a')^3+e{a'-a){a')\^ (2-3) 

where af^^^ = |(ctJ(^) ±2crJ(^)) with , o-^V) > ^|(a) being the usual Pauh spin 
operators at site j corresponding to the a-th {a = 1,2,3) XY spin chain, 
9 {a' — a) is a step function of (a' — a) and ry is a coupling constant; Bi and 
Bl are left and right boundary terms of the form 



3 



2c- exp(2»7) 



exp{r]) 
2c- 



+ sinh T] cosh ^ ^ ^ 
+ sinh^^(J^(i)(J^(2)CT^(3) 



cosh r/ Ea=2 (^t(a) + si^h 7/(7^(2)^1 



(3) 



I 2^^1(3)' 



1 

2C4 



1 



2c^exp(rj) 
I 2c+^L(l)' 



cosh r7Ea=l^i(a) 

+ sinh ?7 cosh ry E^^^ ^ ^ ^ (x£(^) (j£(^) 

+ Smh2 77(7£(i)(7£(2)(^i(3) 

cosh T] ELi f^i(a) + sinh wlii)(^ii2) 



, for m = 1 

for m — 2 
for m = 3 

for 1 = 1 

for 1 = 2 
for 1 = 3 



(2.4) 



(2.5) 



where c± are parameters describing boundary effects. It is worth mention- 
ing that the boundary terms containing three- and two-spin interactions and 
chemical potentials at the left and right ends are consistent with the bulk sym- 
metry which is a combination of the inversion j ^ L — j + 1 and the exchange 
among the sublattices. With the different choices of the pair (m, /) m,l ~ 
1,2,3, it appears that there exist nine classes of integrable boundary terms 
compatible with the intcgrability of the model. After a generalized Jordan- 
Wigner transformation, Hamiltonian (2.1) becomes a strongly correlated elec- 
tronic system with boundary interactions. Quite remarkably, if one restricts 
the Hilbert space to the one which only consists of, say, , cr(2) , then Hamil- 
tonian (2.1) reduces to that of two coupled XY open chains with special 
boundary interactions, which has been considered in [26]. 

We now establish the quantum intcgrability for the system defined by Hamil- 
tonian (2.1), by using the general formalism described in the paper [17]. To 
fix notation, let us briefly recall some basic quantities for the bulk model (2.3) 
with periodic boundary conditions. As was shown in [27], the bulk model 
Hamiltonian commutes with a one-parameter family of bulk transfer matrix 
t{u) of a two-dimensional lattice statistical mechanics model. This transfer 
matrix is the trace of a monodromy matrix T{u), which is defined , as usual, 

by 



T{u) = Lon{u) ■ ■■Loi{u) 



(2.6) 
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with Loj{u) of the form, 



Loj{u) 



-(1), 



(2), 



(2.7) 



where 



a' = l 



+ (^j (a) '^(Ka) + ^tioc) ^O(a) ) 



l + exp(2r) ^ <r+. ,,<r„, ,,)u2. (2.8) 



a'=l 



The commutativity of the bulk transfer matrices t{u) for different values of 
the spectral parameter u follows from the fact that the monodromy matrix 
T{u) satisfies the Yang-Baxter algebra 



12 2 1 

Ri2{u,v) T (u) T (v) (v) T {u)Ri2{u,v) 



(2.9) 



The explicit form of the corresponding R-matrix Ri2{ui,U2) can be found in 
[27]. Here wc only emphasize that the local monodromy matrix as well as the 
quantum R-matrix does not possess the crossing symmetry. It satisfies the 
non-additive Yang-Baxter equation 

Ri2{u, v)Ri3{u, w)R23{v, w) = R23{V, w)Ri3{u, w)Ri2{u, v). (2.10) 

Following the QISM adapted to the case of special boundary conditions, we 
define the doubled monodromy matrix as 

U^u) = T(u)K4u)T-\-u), (2.11) 

such that the boundary transfer matrix is given by 

T{u) = TroK+{u)U_{u), (2.12) 



where T ^ is the inverse of the monodromy T and K± are the matrices defining 
the boundaries. The requirement that they obey the RE [26,28] 



Ruiu.v) {u)R2i{v, -u) K_ {v) 

2 1 

{v)Ri2{u,-v) K_ {u)R2i{-v,-u), 



(2.13) 
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Rl\''{v,u) Kl' {u)R^2{-u,v) K"^ (v) = 

K'^ {v)R2i{-v,u) K'^ {u)R%''{-u, -v) (2.14) 
together with the Yang-Baxter algebra and the following properties 



Ri2{u,v)R2i{v,u) = 1, 
Rl\{-v,u)R%{u,-v) = l, 



-u, v)R 



-u) = l. 



(2.15) 
(2.16) 
(2.17) 



assure that the transfer matrix commutes for different spectral parameters, 
proving the integrability of the model. 

Therefore the transfer matrix (2.12) may be considered as the generating fuc- 
tion of infinitely many integrals of motion for the system. We emphasize that 
there is no isomorphism between the matrices Kj^{u) and K-{u), due to the 
absence of the crossing symmetry for the 7?-matrix. Therefore, we have to 
solve two REs separately in order to fix the boundaries. In the Appendix, 
this calculation is presented in detail using the variable-separation prescrip- 
tion. We find three different classes of boundary K± matrices consistent with 
the integrability of the Hamiltonian (2.1). Let us first fist those related with 
-matrix 



A- 



A_{u) 

B_{u) 

C-(u) 

D^u) 

E_{u) 

F_{u) 

G_{u) 

H^{u) 



,(2.18) 



where for m = 1 



A_(n) = (c_ 


+ M)(e2^'c_ 


+ M)(e'''c_ 


+ u), 


B_{u) = {c_ 


-tx)(e2'?c_ 


+ u)(e^''c_ 




C_(m) = (c_ 


-M)(e2^c_ 


+ M)(e^''c_ 


+ u), 


D_{u) = {c_ 


-M)(e2^c_ 


-M)(e^''c_ 


+ u), 
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E4u) = (c_ - if)(e'^c_ + M)(e^''c_ + u), 
F_{u) = (c_ - M)(e2'?c_ - u){e^'^c_ + u), 
G.{u) = (c_ - M)(e2''c_ - M)(e^"c_ + u), 
H4u) = (c_ - M)(e2''c_ - «)(e^''c_ - u). 



for m = 2 



= E_{u) = (c_ + m)(c_ + e^^'u), 
B_{u)^C4u) = = G4u) = (c_ +ii)(c_ - e^^u), 

D_{u) = = (c_ - m)(c_ - e^^w), 




A.{u) 
B-{u) 



1 



= ^-(m) = (C_ +77,), 



c_ 



However, it is much more tedious to find the boundary K- matrix 
since not only the corresponding RE is more involved but also the new object 
i?i2(M, is more complicated. We list the final result here and some details 
can be found in Appendix, 



AM 


















B+{u) 












C+{u) 












D+{u) 












E+{u) 












F+{u) 












G+{u) 












H+{u) 



for / = 1 
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A+{u) = (e^^c+u - l){e^''c+u - l){e''''c+u - 1), 

B+{u) = e^'iie^'^c+u + l){e^'^c+u - l){e^'^c+u - 1), 
C+{u) = e^'\e^'hi^'u + l)(e'^^'c+i/, - \){(?''^c^u - 1), 
D^{u) = e^\e^"c+u + l){c+u + l){e^'^c+u - 1), 
E+{u) = {e^"c+u + l)(e^''c+M - l){e^''c+u - 1), 
F+(7i) = e^"{c+u + l)(e2''c+n + l){e^"c+u - 1), 
= e^"{c+u + l)(e2''c+n + l)(e2''c+w - 1), 
H+{u) = e^''(c+« + e^'^){c+u + I)(e2'?c+M + 1). 

for / = 2 



A+{u) = = {e^^'c+u - l)(e^^c+M - 1), 

C+{u) = = e^'^ie^^'c+u + l)(e^^c+M - 1), 

E+{u) = = (e^'^c+M + l)(e^''c+M - 1), 

G'+(w) = H+{u) = e2''(c+w + l){e^''c+u + 1), 

for / = 3 



^+1 




= B+ 


(u) 


= e'^'(e^''c+n - 1) 




[^) 


= D^ 


.{u) 


= (e^V«-l), 




[u) 


= F+ 


(u) 


= e2''(c+« + l), 




[u) 


= H^ 


-H 


= (c+w + 1), 



The above explicit formulae for K±{u), derived by solving the two REs directly, 
clearly show that no automorphism between K^{u) and K^{u) exists and 
K^{u) can not be obtained from K_{u). These three classes of boundary K±- 
matrices provide nine possible choices of BC, according to the combination of 
the boundary pairs (K^\u), K^\u)), m,l — 1, 2, 3, which originate the pair- 
boundary terms (Bi^\ B^^)^ respectively. Taking into account the fact that 

TrKX'i^) ^0 ,Tr K+ (0) = 0, Tr K+ (0) = and TrK^+>{0) = 0, Tr K+ 
(0) = 0, as well as TrK^^\o) = 0, we can show that the Hamiltonian (2.1) 
with the boundary terms pairs {K^\u) , K^\u)) is related to the transfer 
matrix matrix (2.12) in the following way 



t{u) = Ciu^ + C2{H + const.) u'^ + ■ 
t{u) = C^u^ + Ca{H + const.) + • 
t{u) = Cf,u + C(y{H + const.) it^ + • • 



for / = 1, (2.20) 
for I = 2, (2.21) 
for I = 3. (2.22) 



Above Q, i = 1, • ■ - 6, are some scalar functions of the boundary parameters. 
The symbols, either "bullet" or "prime" denote the derivative with respect 
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to the spectral parameter. For 1 = 1, the Hamiltonian (2.1) is related to the 
fourth derivative of the boundary transfer matrix t{u) which can be derived 
by the simplified formula [28] 



8TrK^+> (0) j=i ^ 

^ TT{Kf{0)L',^{0)PoL), (2.23) 



Tr^i'^ (0) 

Hjj^, = Po,J^^moMPo'J\0)P,-j^M (2.24) 
with m = 1, 2, 3. 

However, to obtain new boundary terms corresponding to the X^'^ -matrices 
when I = 2,3, we have to develop new formulae regarding to the relations 
(2.21) and (2.22). What follows are the general formulae of the Hamiltonian 
related to Eqs (2.21) and (2.22) given by 



+ j {lV(i^f (O)Lol(O)Po^) +TV(i^f (O)PolZoI(O)) 

+6 IV(i^f "(O)Lo^(O)Pol) + 3 Tr(Kf '(O)C(O)Pol) 
+3 IV(xf '(0)Po^l;'JO)) + 3 'IV(Kf (0)l;JO)ZoJO)) 
+3 TV(i^f (0)LoJO)L;JO)) + 6 Tr(Kf '(0)L;^(0)4J0)) \ , (2.25) 



and 



L-l 



= E ^..+1 + ^^oii^i"^^ (0)Pc 



Tri^rl*^ (0) i=i 



+ 



{Tr(Kf (O)C(O)Pol) + Tr(Kf (0)Pol4'l(0)) 



Tri^i^^ (0) 

+2 Tr(irf (0)4^0)4^0)) + 2 Tr(Kf '(0)LoJO)Poi) ^ , (2.26) 



respectively. Above 



4l(0) = ^Loi(-M)U=o = PoLi^'(0)PoL, 
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4'l(0) ^ ^,L,l{-u)U=o + PolL"{0)Pol, 



/ = 6 



Tr(i^f (0)LoJO)Pol) + Tr(i^f (0)P^o4'l(0)) + TrK^^ (0|l27) 



Using above formulae and after lengthy and tough calculation, one can derive 
the boundary terms presented by Eqs (2.4) and (2.5). It was seen that the 
Hamiltonian of the model can be derived from expansion of the boundary 
transfer matrix (2.12) at m = up to the fourth, third and second orders 
respectively. This contrasts to the case of two coupled XY spin chain where 
the Hamiltonian with four classes of the boundary terms coming from the third 
and second derivatives of the transfer matrix [26]. In the next sections, we shall 
be focusing on the solution of the eigenvalue problem of the Hamiltonian (2.1). 



3 The Bethe ansatz equations 

Having established the quantum integrability of the model, let us now solve it 
by using the coordinate space Bethe ansatz method. Following [19], we assume 
that the eigenfunction of the Hamiltonian (2.1) takes the form 

|V&)= ^ ^'a„...,a.<„r--<a.|0)' (3-1) 

{{xj,aj)} 

N 

*ai,-,a;v(^l> • • • > ^iv) = (^pAaQi,.-,aQN{kpQi, kpQN) exp(i ^ kp.Xj), 

P j=l 

where the summation is taken over all permutations and negations of A;i, • • • , k^, 
and Q is the permutation of the particles such that 1 < xqi < • • • < 
xqn < L. The symbol ep is a sign factor ±1 and changes its sign under 
each 'mutation'. Substituting the wavefunction into the eigenvalue equation 
H\<if) = one gets 

-^■■■,aj,ai,--{' ' ' ) kj, ki, • • ') — Sij{ki, kj)A...^a^^o,j,— {' ' ' 1 kit kj, • • •), 

-Aai,-{ — kj, ■ ■ •) — S^{kj■,Pla^)Aa^^...{kj, ■ ■ •), 

A..,a,(- • • , -kj) = s''(kj;pLa,)A...,c.,(- ■ ■ , kj) (3.2) 
with Sij{ki, kj) = Sij{^{ki — kj)) being the two-particle scattering matrix, 

sin k 
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5-(*,^-ie<~-;^, a^P, (3.3, 
and s'"{kj\piai) ^^i^j'tPLai) the boundary scattering matrices, 



s'^ik.^p.J = iZ^^-^^e^«^^(^+^) (3.4) 
1 - Pr.rv.e*'^^ 



where Ofj = 1,2,3 and pi^j and p^ai are given by the following formulae, 
corresponding to the nine cases with respect to the boundary pair {b[™'^ , Bj^^), 
respectively. 



Case i : pu = pu = Pis = Pi 



c_ exp{4r)) ' 



m=P.2=P.3-PL = ^^^^; (3.5) 

Case ii : pu = 0, pu = Pia = Pi+ = — , 

c_ 

Pli = Pl2 = Pl+ = ^ exp(2?7) ' ^ °' ^^'^'^ 

Case iii : pu = P12 = 0, P13 = —, 

c_ 

Pli = — , Pl2 = Pl3 = 0; (3.7) 

Case iv : pn = pi2 = P13 = Pi = 

c_ exp(4?7j 

1 , , 

Pli = Pl2 = 777^, Pls = 0; (3.8) 

c+exp(2?7) 

Case V : pn = pi2 = P13 = Pi 



c_ exp(47y) ' 

1 , , 

Pli = — , Pl2 = Pl3 = 0; (3.9) 

Case vi : pn = 0, pi2 = P13 = — , 

c_ 

1 

Pli = — , Pl2 = Pl3 = 0; (3.10) 
c+ 

Case vii : pn = 0, pi2 = P13 = — , 

c_ 

Pli=Pl2 = Pl3^Pl = ^^^^; (3.11) 
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Case viii : 



PU = Pl2 = 0, pi3 = — , 



PLI = PL2 = PL3 = PL 



1 



c+ exp(2?7) 



(3.12) 



Case ix : pn = = 0, pig = — , 



1 



Pli = PL2 = 



c+ exp(2?7) 



, Pl3 = 0; 



(3.13) 



As is seen above, the two-particle S-matrix (3.4) is nothing but the R matrix of 
the A2 XXZ Heisenberg model (in the homogeneous gauge) and thus satisfies 
the quantum Yang-Baxter equation, and the boundary scattering matrices 
and obey the corresponding reflection equations (for details, see, [29]). This 
is seen as follows. One introduces the notation 

s{hp) ■ (3.14) 

1 — pe^^ 



Then the boundary scattering matrices s^{kj;piai), s^{kj;pLai) ^e writ- 
ten as, corresponding to the nine cases, respectively. 



Case i : s^{kj;piaj = s{-kj;pi)I, 
s''{kf,pLa,)^e^'^'^'^+'^s{k,;pL)I; 

Case ii : s^{kj;piai) = s{-kj;pi+) 



(. 



s^{kf,pr.a.)^e''^'^^+'h{k,;pr.,) 



Q ^iki sin(iC-+fcj/2) 
sin(j(;^- — fcj/2) 

1 
1 

Q Q j-.iki sin(m++fcj/2) 

sin(iK+— fcj/2) / 








(3.15) 



iki sin«-+fcj/2) 
sin(iC--fej/2) / 



(3.16) 



Case iii : s^{kj;piai) = ^(-^jSPii) 



1 
1 






jfe . sin(m_+fcj/2) 



e^'^J ^-TT 



(m_-fcj/2) / 



«^(%;PLaJ = e''=^'^(^+^).(A;,;p^i) 



1 







Q ^ifcj sm(iC++fcj/2) 
sin(if+— fej/2) 







1(9.17) 

gifej sin«++fcj/2) 
sin(if+— fcj/2) / 
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Case iv : s^{kj]piai) = s{-kj;pi)I, 

( 



Case V : s^{kj;piai) = s{-kj]pi)I, 

(. 



1 

1 









(3.18) 



QQ^,fc^. sin(«++fc,/2) 

sin(i/t+— fej/2) / 












sin(if_|.— fcj/2) 
Q Q gifcj sin(iC++fcj/2) 



Case vi : s^{kj\pxa^ = s{-kj\px+) 



(. 



s''{kj;pL^^)^e^'^-'^^^'^s{kf,pLi) 



,ikj2{L+l) , 





g ^ikj sm(tC-+fcj/2) 
sin(j(^_ —kj/2) 






sin(if_l_— fej/2) y 





(9.19) 





ikj sin«-+fcj/2) 
sin(iC--fej72) / 



g gtfcj sin(iC++fcj/2) 
sin(if+— fej/2) 










^jfcj sin(iC++fcj/2) 
sin(i(j4-— fcj/2) / 



Case vii : = s{-kj;pi+) 

s''{k,;pLa,) = e'''^'^'^+'h{k,;pL)I; 
Case viii : s^{kj;piai) — "^("^jiPii) 
«''(A^,;PLaJ = e^'=^-^(^+^).(%;p^)/; 
Case ix : s^{kj;piai) = 



1 







g ^ikj sin(iC-+fcj/2) 
sin(if_— A;j/2) 










(9.20) 



ikj sin(iC-+fcj/2) 
^ sm(iC--fcj/2) / 



(3.21) 



1 
1 








g g ^iki sin(m-+fcj/2) 

sin(m_— fej/2) / 



(3.22) 



1 
1 








g g ^ikj sm(m_+fcj/2) 



1 

1 



sin(iK;_— fej/2) / 

^ 





g g gifcj sin(iK++fcj/2) 



(3.23) 



sin(iK;_|_— fej/2) y 

Here / stands for 3 x 3 identity matrix and pl+ are the ones given in 
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(3.6); (^±, K± are parameters defined by 



(3.24) 



We immediately sec that (3.15) are the trivial solutions of the reflection equa- 
tions, whereas (3.16) and (3.17) are the diagonal solutions for the A2 XX Z 
model R matrix. The boundary scattering matrices given in (3.15), (3.16) and 
(3.17) also constitute nine classes of possible choices of boundary conditions 
for the spin degrees of freedom of the model. 

Then, the diagonalization of Hamiltonian (2.1) reduces to solving the following 
matrix eigenvalue equation 

Tjt = t, j = l,---,N, (3.25) 

where t denotes an eigenvector on the space of the spin variables and Tj takes 
the form 

Tj = 57(fc,)s^(-A;,;pi„j7?7(/c,)i?+(/c,)5^(A;,;pi«J5+(A;,) (3.26) 

with 

~ ^j,N{kj, k]\i) ■ ■ ■ Sjj^i{kj, kj^i), 
{kj) = Sjj-i{kj, kj-i) ■ ■ ■ Sj^i{kj, ki), 
Rj {kj) = Sij{ki, —kj) ■ ■ ■ Sj-ij{kj-i, —kj), 

R^ikj) = Sj+ij{kj+i, -kj) ■ ■ ■ SN,j{kN, -kj). (3.27) 

This problem may be solved using the algebraic Bethe ansatz method. The 
Bethe ansatz equations are 



e^'^^'(^^')F{k,;pu,PLi)=l[ 



sin[i(A:,-AW) + f]sin[i(fc,+AW) + f] 



a=i sin[i(A;,- - aL')) - |] sin[i(A;,- + AL' 



(l)^ 



N 

n 



sin[i(A« - k^) + f ] sin[i(A« + A:„ + f ] 



isin[i(A?)-A;„)-f] sin[i(A 



(1) 



121 
2 J 



X 



M2 

n 

5=1 sm[f 



sin[i(A(^) 



ir7]sin[i(A?^ + A:;: 



■7' ) 



sin[i(AW-Af)-f]sin[i(AW + Af)- 



) + f ] 



+ ^|sin[i(A 



(1) 

7 



+ ir]] 



it]] 
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sm[i(Ap - A^^^) + ///] sm[i(Af + A^^^) + ii]] 



n sin[|(A?) - A^)) - irj] sin[|(A?) + A^)) - irj] 
sin[i(Af - AW) + f ] sin[i(Af + A^) + f ] 

sin[i(Af - Ai^)) - f ] sin[|(Af + A^) - f ] ' ^ ^ 

where 



F{kj;pi+,pL+) = s{kj;pi+)s{kj;pL+), (for all cases) 



G'(AW;C-,C+) = < 



i^(A(2);K_,«^) = < 



sin(-iC- -iAW -7/2+rj/2) tA^^^ 
sm(iC-+iA(^' -7/2+57/2) 



sin«+-iA(i) -7/2+77/2) 
sin(iC++iA(i) -7/2+77/2) ' 

1 

siii( 



(1) 



sin(' 
sin(' 



sin(' 

sin( 



:^^_iA(i)_^/2+,,/2) iA' 
;C++jA(i)-7/2+T;/2)^ 
:C- -tA(i) -7/2+77/2) siii«+-iA(i) -7/2+77/2) 



:C-+iA(l)-7/2+r;/2) sin(7C++iA(l) -7/2+77/2) 
g_-iA(i)-^/2+77/2) a'i) 



sin(7C_ +iA(i) -7/2+77/2) ^ 
1 



1 

sin(jK-|-— I. 



A(: 



sin(i«;++iAP 
sin(iK_ — iA(^ 



sin(iK_+iA(- 
sin(iK+— iA(^ 



sin(iK_|_ +iA(2^ 
1 

1 



-7/2+77) ^-»A<,'^ 



-7/2+77) 
17/2+7?) giA^2) 



-7/2+77) 
-7/2+T?) ^-iA^ 



(2) 



-7/2+77) 



case 
case 
case 
case 
case 
case 



case (i) 
case (ii) 
case (Hi) 
case {iv) 
case (v) 
case (vi) 
case (vii) 
case {via) 
case (ix) 

i) 
ii) 

Hi) 

iv) 

v) 

vi) 



1 

sin(iK- 



-zAt'^) -7/2+77) ^iA 



(2) 



sin(7K_+7A(2) -7/2+77) 
sin(iK_ — 7A(-^) —7/2+77) sin(m+— jA(^) -7/2+77) 
sin(iK;_+iA(2)— ')//2+T7) sin(iK;++iA(2)— 7/2+77) 



case 



case 



case 



Vll) 

via) 

ix) 



(3.29) 



The energy eigenvalue E of the model is given by = — 2 Y^^=i cos kj (modular 
an unimportant additive constant). 
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4 Conclusion 



In this paper, we have studied integrable open-boundary conditions for the 
three coupled XY spin chain modeL The quantum integrabihty of the bound- 
ary system has been estabhshed by the fact that the corresponding Hamil- 
tonian may be embedded into a one-paramctcr family of commuting transfer 
matrices. A desirable way to handle open boundary conditions for the mod- 
els associated with a general class of quantum R- matrices (with or without 
crossing-unitarity) has been developed. Moreover, the Bethe Ansatz equations 
are derived by means of the coordinate Bethe ansatz approach. This provides 
a basis for computing the finite size corrections to the low-lying energies in the 
system, which in turn could be used together with the boundary conformal 
field theory technique to study the critical properties of the boundaries. 

Lastly, it is interesting to formulate a graded version of the quantum Yang- 
Baxter algebra and reflection equation algebra for the fermionic Bariev model. 

The algebraic Bethe ansatz for the three coupled XY model with both periodic 
and open boundary conditions would be very significant in understanding of 
the symmetry structure of the model. Those will be addressed in a future 
publication. 
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A Derivation of the boundciry K-t-matrices 



The R- matrix is a 64 x 64 matrix with 216 non-zero elements which are given 
in [27]. We now are looking for diagonal solutions K±{u) of the REs. We 
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parametrize K±{u) as 



zl±{u) 
z2±{u) 
























z±3{u) 












zA±{u) 













z5±{u) 













z6±{u) 












z7±{u) 




^ 









zS±{u) J 



.(A.l) 



Using some simpler functional equations from the first RE wich are given by 
(A.1)-(A.8) in the Appendix A of the paper [28], one can present an ansatz 



z2_ 


[u) 


_ Ci - c^u 


z3_{u) 


C2 - CgU 


(A.2) 


zl_ 


[u) 


Ci + CgU ' 


zl_{u) 


('2 + 1^9 ' 


2:5 _ 


'u) 


_ C3 - CloU 


zA-{u) 


_ C4 - CiiM 


(A.3) 




u) 


C3 + CloU 


z2_{u) 


C4 + Ciiu' 




u) 


C5 — C12U 


z7-{u) 


C6 - Ci3li 


(A.4) 


z2_ 


u) 


C5 + C12U ' 


z3-{u) 


C6 + C13M ' 


z8_ 


u) 


C7 — C14U 






(A.5) 


z4_ 


u) 


C7 + CuU 







with minimal coefficients q to be determined. Running the RE again with 
above ansatz, it can be found that only one coefficient is free and three classes 
of boundary i^_-matrices can be immediately chosen as the forms presented 
in (2.18). 

To solve the second RE is rather cumbersome and sophisticated. For clarifying 
the functional equations arising from the second RE, we denote the Boltzmann 

weights associated with the i?-matrix as W\[u^v\ - • ■ ^ Wr^^{i.L, v) in accordance 
with the orders listed in the paper [27]. The convenient notations z = z^{v) 
and z = z+{u) will be implied hereafter. Similarly, after substituting K^- 
matrix (A.l) into the RE (2.14), we may pick up some simpler functional 
equations such as 



z2 u!2{v, u)pi{v, —u)z2 + w^{v, u)p2{v, —u)zl 

zl w^{—u, —v)p2{u, —v)e~'^'^z2 + W2{—u, —v)pi{u, —v)zl ' 



(A.6) 
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z3 W2{v, u)pi{v, —u)e'^^z3 + w^i^v, u)p2{v, —u)zl 

zl W3{—u, —v)p2{u, —v)z3 + W2{—u, —v)pi{u, —v)e^'^zl ' 

zb u]2{v, ii)pi{v, —u)z5 + w^{v, u)p2{v, —u)e~^^zl 

zl W3{-U, -V)p2{u, -V)zb + W2{-U, -v)pi{u, -v)zl ' 

58 W4:j{—u, —v)p3{u, —v)e^^z6 + W4q{—u, —v)p^{u, —v)zS 
56 W4Q{v,u)pi{v,—u)zQ + W4:j{v,u)p2,{v,—u)e'^'^z'& ' 

56 wi^{—u, —v)pq{u, —v)e^^z2 + W2o{—u, —v)p^{u, —v)e^'"^'^zQ 
52 w\Q{v,u)pQ{v,—u)e'^^zQ-\-W2Q{v,u)p^{v,—u)harz2 

54 _ W2o{—u, —v)p5{u, —v)z2 + wig{—u, —v)p'j{u, —v)zA 

52 wig{v, u)p'j{v, —u)z2 + W2o{v, u)pz{v, —u)zA 

57 Wig{—u, —v)py{u, —v)e^^zl + W2q{—u, —v)p5{u, —v)z3 

53 Wig{v, u)p'j{v^ —u)e^'^z3 + W2q{v, u)p5{v, —u)e^^ harzl 



(A.7) 
(A.8) 
(A.9) 
(A.IO) 
(A.ll) 
(A.12) 



Before further going on, we first need to work out the matrices .^21*^ (^j ^) 1 
Ri2^{u,v), Ri2{u,v) and R2i{u,v) according to their definitions in section 2. 
For our convenience, we prefer to present below some entries of R12 involving 
in above equations as 



{l + e^''uv){l + e^'iuvf{l + e^'^uv){l + uv) 



(1 + e^''uv){l + e^''uv){uv + e^'^){l + uvf^{l + u^){l+v^) 

^4^^ =^7 T^A;^ TTY^ ^-^^ 



{l + e^iuv)\l + e^''uv){l + uv)^ 

Pn^^'")- (^^_^y(^^2vu-v){e^'^v-u){e'^'^v-u)e'^' ^ ' ^ 

Then by analying the structure of the Eqs. (A.6)-(A.12), we can get the fol- 
lowing relations: 



z2^{u) ciu + e^cs 
zl+{u) 
z5^{u) 



e?^C2U — Cg ' 

zl^{u) e^'^c^u — Cio' z2^{u) c^'^CqU — Cn ' 



CiU - Cs 
C3U + Cio 



z3+{u) + e^'^Cg 

zl+{u) 
zAju) 



(A.20) 
(A.21) 
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z6+{u) ^ C5U + C12 z7+{u) ^ C7U + Ci3 

^6+(m) e'^^CiU — Ci4 

Running second RE (2.14) with above relations again and again, the solutions 
(2.19) would be fixed definitely. 
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